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Abstract 

We have investigated the pressure induced structural and electronic properties of plutonium 
pnictides (PuY, Y= P, As, Sb). The total energy as a function of volume is obtained by means of 
self-consistent tight binding linear muffin-in-orbital (TB-LMTO) method within the local spin 
density approximation (LSDA). From present study with the help of total energy calculations 
(spin polarized) it is found that PuP, PuAs and PuSb are stable in NaCl - type structure under 
ambient pressure. The structural stability of PuP, PuAs and PuSb changes under the application 
of pressure. We predict a structural phase transition from NaCl-type (Bi-phase) to CsCl-type 
(B 2 -phase) structure for these Pu-pnictides in the pressure range of 20.8 - 42.0 GPa. We also 
calculate the lattice parameter, bulk modulus, band structure and density of states. From energy 
band diagram it is observed that all the three compounds exhibit metallic behaviour. The 
calculated equilibrium lattice parameters and bulk modulus are in good agreement with available 
experimental data. 
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1. Introduction 

The rare earth lanthanide and actinide binary compounds exhibit variety of solid state properties 
due to presence of partially filled /-electrons in these systems. The nature of /-electrons, 
however, can be changed from localized to itinerant, leading to a variety of changes in their 
structural, magnetic and electronic properties under the application of pressure [1-2]. The 
pressure induced structural phase transition in binary compounds of lanthanides and actinides 
with NaCl-type structure has received considerable attention in the recent past, due to their 
extensive applications in different fields [3-4]. In the case of actinide compounds, the RE ion has 
electronic configuration 5f (6d7s) m . Under the application of the pressure, these compounds lead 
to structural changes associated with a valance change to 5f' 1 {6d7s) mJtl . Since cohesion of solids 
is crucial, the interionic distance changes as a consequence of the application of pressure. 
Therefore the derealization or itinerant behaviour strongly depends on actinide-actinide distance 
under pressure. Recent studies using synchrotron and powder X-ray diffraction demonstrated that 
under pressure, a majority of binary compounds with a NaCl-type (Bi) structure undergo a phase 
transition to CsCl-type (B 2 ) structure with eightfold co-ordination at high pressure [5-20]. 
Considerable theoretical but few experimental studies on plutonium binary compounds have 
been reported during the last few years. In the mono-pnictides of plutonium PuAs shows Bi to B 2 
transition around 35 GPa [1] while PuSb undergoes a B 2 transition at 18 GPa [6]. The structural 
and elastic properties of PuAs and PuSb compounds have been theoretically studied by 
Srivastava et al [5], using an interionic potential model approach. No such effort has so far been 
made in the case of Plutonium phospide. Moreover, spin polarised and pressure induced 
electronic structure and its correlation with the structural phase transition in is yet to be explored. 

It is well known that with increasing atomic number, the /-shells become smaller and thus 
are more localized implying a smaller band width [7]. This localization is also favoured from 
lighter (e.g. P) to heavier anions (e.g. Bi). In this context, it will be interesting as well to 
investigate spin polarized structural phase transition and electronic properties in PuP and 
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correlate the same with electronic properties for a few selected Plutonium Pnictides, with fixed 
atomic number but changing the anion mass (P to Sb). 

In the present paper we have investigated, for the first time, the high pressure structural and 
electronic properties of PuP at ambient and high pressure using ab initio Tight Binding-LMTO 
method within the local spin density approximation (LSDA). Besides, we also report similar 
results for PuAs and PuSb and compare their structural and electronic properties to visualize how 
the electronic properties are related to each other in this particular group. In the section 2, we 
briefly describe the method of total energy calculation which will be followed by results and 
discussion on structural and electronic properties in section 3. 

2. Method of Calculation 

The total energy, structural stability, electronic band structure (BS) and density of states (DOS) 
of three Plutonium pnictides (PuY, Y =P, As and Sb) are calculated in the non-magnetic (NM) 
and ferromagnetic (FM) states by using the tight binding LMTO (TB-LMTO) method [24-25] 
within local spin density approximation (LSDA) [26]. Von - Barth and Hedin parameterization 
scheme [27] has been used for exchange correlation potential, similar to our previous work [23]. 
The density of states was calculated by tetrahedron method [28]. All the PuY compounds 
considered here crystallize in Bi-type (NaCl) structure (Space group Fm 3 Dm, No 225) and 
stable in ferromagnetic state at ambient pressure. In Bi structure, Pu atoms occupy (0,0,0) and Y 
atoms at (0.5, 0.5, 0.5) positions respectively. Since TB-LMTO method works fairly well for 
closed packed structures and since PuY compounds in NaCl structure are not closed packed, 
therefore, two equivalent empty spheres were introduced at positions (0.25, 0.25, 0.25) and 
(0.75, 0.75, 0.75), in such a way that they do not destroy the crystal symmetry [29]. However in 
the case of CsCl- type structure at high pressure, it being a closed packed one; no such empty 
spheres were needed. The calculations were performed for 8000 k- points (20x20x20 k 3 grids) in 
the Brillouin zone for both Bi and B 2 structures. The E and k convergence were checked 
subsequently. The total energies were calculated by changing the volume from 1.05 to 0.65 V , 
where V is the equilibrium cell volume at ambient pressure. The calculated total energies as a 
function of volume for AnP compounds in Bi and B 2 structures were fitted to Birch equation of 
state [30] to obtain pressure - volume(P-V) and pressure-enthalpy (P-H) relations, and finally 
leading to equilibrium lattice constants (a) and bulk modulus (B) . 
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3. Result and discussion 

3.1. Structural properties 

The magnetic stability of PuY compounds has been investigated by estimating the total 
energy as a function of volume in NaCl structure by performing spin-polarized (FM) and non - 
spin polarized (NM) calculations. Fig.l (a-c) shows the calculated total energy of these 
compounds in FM and NM states. It is revealed from this figure that PuP, PuAs and PuSb are 
stable in FM state with minimum energy.The structural stability of PuY (Y = As, Sb and P) 
compounds has been examined by calculating the total energy as a function of volume in 
ambient (NaCl phase) and at the high pressure (CsCl phase) structures in FM state is plotted in 
Fig.2 (a-c). Fig 2. clearly indicates that these compounds are stable in Bi (NaCl) phase at 
ambient and there is a possibility of structural phase transition to B 2 (CsCl) phase at sufficiently 
reduced volume (or high pressure). These compounds, as is seen from Fig. 2, transform to CsCl 
structure on compression. The calculated values of equilibrium lattice constants, in the NaCl 
phase are 5.51, 5.76, 6.18 A for PuP, PuAs and PuSb, respectively which agree fairly well with 
the experimental data [7], shown in Table 1. The structural phase transition is determined by 
calculating Gibb's free energy (G) for two phases (Bi and B 2 ) which is given by G=E to t + PV - 
TS. Since the theoretical calculations are performed at T= OK, Gibbs energy becomes equal to 
enthalpy (H) (H = E tot + PV). For a given pressure, a stable structure is one for which enthalpy 
has lowest value and the transition pressure (P t ) is determined at which the enthalpies of two 
structures are the same. Following such a procedure, the phase transition pressures P t for PuY 
compounds are calculated for a Bi - B 2 transition. The phase transition pressures have been 
estimated as 42 GPa, 37 GPa and 21 GPa, for PuP, PuAs and PuSb with volume collapse of 
8.9%, 9.4% and 8.6% respectively. The pressure - volume relationship or equation of states for 
these compounds are shown in Fig. 3 (a-c). The insets show the variation enthalpy as a function 
of pressure for respective compounds. In the B 2 phase the corresponding lattice parameters are 
3.32, 3.51 and 3.71 A respectively. In addition to the above structural properties, the values of 
bulk modulus (B) in both the structures are calculated by using Birch equation of state [30], 
which has been shown in Fig.3 (for energy as a function of pressure). The calculated results are 
presented in Table 1 and compared with available experimental results for all the three 
compounds [5,7,19]. Our calculated values are, however little larger, as compared to the 
experimental values, the reason being the same as in the case of lattice parameter. However a 
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more precise measurement might be helpful looking to the uncertainties in the experiment data. 
No such experimental results in the B 2 phase are reported so far and could not be compared at 
present. 



3.2 Electronic Properties 

To understand the electronic properties, we have carried out spin - polarized electronic band 
structure (BS) and DOS calculations of PuY compounds in both Biand B 2 phases. As mentioned 
earlier, these compounds are metallic ferromagnet. We therefore, have plotted the band 
structures and density of states in majority and minority spin channels at ambient (Bi phase) as 
well as at high pressure (B 2 phase). An exact understanding of the electronic structure of PuY 
compounds is an extremely challenging problem, because of the itinerant behaviour of their 5-/ 
electrons. However a qualitative picture of their electronic properties can be drawn by comparing 
the same different pnictides of Pu-element. We have, therefore, calculated the electronic band 
structure (BS) and density of states (DOS) of three PuY compounds in their Bi and B 2 phases 
and plotted them in Figs. 4(a-c) and 5(a-c), where the solid line represents the BS and DOS of 
majority spin while dotted line represents BS and DOS for minority spin at ambient pressure. 
The Fermi level is shown by single horizontal dotted line and is set to zero. It is seen from Fig 4 
that the PuY compounds are metallic in nature in Bi phase at ambient, and the region around 
Fermi energy is densely populated by Pu-d and /- bands .To understand the band structure of 
PuY qualitatively, in Figs 6 to 8, we have plotted the partial density of states (P-DOS) . These 
figures reveal that lowest bands around -1.4Ry. are due to Pu-s bands which hybridize with Y-s 
bands in all the three pnictides . The Pu-d band is quite broad and spread over the energy 
interval -0.75 to 1.0 Ry. highly localized Pu-/band, stays above the Fermi level, at ambient and 
moves marginally towards Fermi level on application of pressure (not shown in P-DOS for B2 
phase). 
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4. Conclusion 

The structural phase transition due to high pressure and electronic properties of PuY compounds 
have been studied using TB-LMTO method. All the three compounds undergo a structural phase 
transition to B2 structure at 37, 29 and 24 GPa, respectively. All the compounds are found to be 
metallic in both the phases. 

Acknowledgements 

The authors are thankful to Dr. V Srivastava, for useful discussions and CSIR, and UGC (SAP) 
for financial assistance. 

References 

[I] S. Dabos-Seignon, U. Benedict, S. Heathman, J.-C. Spirlet, M. Pages, J. Less Common Met. 160 
(1990) 35-52. 

[2] L. Petit. A. Savane, W. H. Temmerman, Z. Szotek, Euro. Phys. J. B 25 (2002) 139-146. 

[3] Chun-.Gang. Duan, R.F. Sabirianov, W.N. Mei, P.A. Dowhen, S.S. Jaswal and E.Y. Tsymball, 

J. Phys. Cond. Matter 19 (2007) 315220 -315251. 
[4] CM. Aerts, P. Strange, M. Home, W.M. Temmerman, Z. Szotek and A. Svane , 

Phys. Rev. B 69 (2004) 0451 15 -045120. 
[5] V.Shrivastava, S. P. Sanyal, J. Alloys Comp. 366 (2004) 15-20. 
[ 6] S. Dabos-Seignon, U. Benedict, High Pressure Res. 4 (1990) 384-391 

[6] P. Burlet, J.M. Fournier, E. Pleska, S. Quezel, J. Rossat-Mignod, J.C. Spirlet, J. Rebizant and O. 

Vogt, Journal De Physique Colloque C8, 12 (1988) 469-474. 
[7] O. Vogt and K. Mattenberger, J. Alloys and Comp. 223 (1995) 226-236. 
[8] A. Singh, V. Srivastava, M.Aynyas and S.P. Sanyal, J. Nucl. Mater. 401 (2010) 60-64. 
[9] P. M. Oppeneer, T. Kraft , M.S.S. Brooks , Phys. Rev. B 61 (2000) 12825-12830. 
[10] J. S. Olsen, L. Gerward, U. Benedict and S. Dabos, O. Vogt, Phys Rev B37 (1988-11) 8713-8718. 

[I I] L. Gerward, J. Staun Olsen, U. Benedict, S. Dabos, J.-P. Itie, O. Vogt, High Temperatures High 
Pressures 20 (1988), 545-552. 

[12] L. Gerward, J. Staun Olsen, U. Benedict, S. Dabos, O. Vogt, High Pressure Res.l (1989) 235- 
251. 



A Quarterly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories 
Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., [•MJiPfleMJFJILMM as well as in Cabell's Directories of Publishing Opportunities, U.S.A. 

International Journal of Engineering, Science and Mathematics 
http://www.ijmra.us 




ISSN: 2320-0294 



[13] M. Gensini, E. Gering, S. Heathman, U. Benedict, J.C Spirlet, High Pressure Res.2 (1990) 347- 
359. 

[14] P.K. Jha, S.P. Sanyal, Phys. Rev. B 46 (1992) 3664-3667. 

[15] Fournier, J.M. Troc, R., Handbook on the physics and chemistry of actinides, Eds A.J. Freeman 

and G.H. Lander (North Holland Amsterdam) 2 (1985) , 29. 

[16] V. Srivastava, M. Rajagopalan S. P. Sanyal, J. Mag. Mag. Mat. 321 (2009) 607-612. 

[17] P. Pandit, V. Srivastava, M. Rajagopalan, S. P. Sanyal, Physica B 403(2008) 4333-4337. 

[18] U. Benedict, J. Less-Common Metals 128 (1987) 7-45 

[19] U. Benedict, /. Alloys Comp. 223 (1995) 216-225 

[20] Y. Meresse, S.Heathman, C. Rijkeboar, J. Rebizant, J.Alloys Comp. 284 (1990) 65-69. 

[21] B.S. Arya, M. Aynyas, S.P. Sanyal, J. Nucl. Mater. 393 (2009) 381-386. 

[22] J.H. Handwerk , O.L. Kruger, J. Nucl. Engineering and Design, 17 (1971) 397-408. 

[23] Chandrabhan Makode and Sankar P. Sanyal, Physica B (201 1) 3 175-3 179. 

[24] O.K. Andersen, Phys. Rev. B 12 (1975) 3060-3083. 

[25] O.K. Andersen, O. Jepsen, Phys. Rev. Lett. 53 (1984) 2571-2574. 

[26] W. Kohn, L.J. Sham, Phys. Rev. A 140 (1965) 1 133-1 138. 

[27] U. van Barth, L. Hedin, Phys. C 5 (1972) 1629-1642. 

[28] O. Jepson, O.K. Andersen, Solid State Commun 9 (1971) 1763-1767. 

[29] N.E. Christensen, Phys. Rev. B 32 (1985) 207-228. 

[30] F. Birch, J. Geophysics, Rev. 83 (1978) 1257-1268. 



A Quarterly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories 
Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., [•MJiPfleMJFJILMM as well as in Cabell's Directories of Publishing Opportunities, U.S.A. 

International Journal of Engineering, Science and Mathematics 
http://www.ijmra.us 




ISSN: 2320-0294 




Table 1: Calculated equilibrium lattice parameter ao (A), bulk modulus B (GPa,), phase transition pressure 
P t (GPa), Magnetic moments juB and relative change in volume at transition pressure of PuP, 
PuAs and PuSb in Bi and B 2 Phases. 



Solids Lattice Bulk Modulus Magnetic Phase Relative 

Constant B (GPa) moment (juB) Transition Volume 

o 

(A) Bi B2 Bi B2 Pressure Change 

Bi B 2 (GPa) (atPt)% 



PuP 

Present 5.51 3.32 146.2 112.6 4.8 4.9 

Experimental 5.55 a — 

Other 

PuAs 

Present 5.767 3.51 63.9 56.2 4.5 5.01 

Experimental 5.88 a — — — — 

Other 5.84 b — 69 b 

PuSb 

Present 6.187 3.75 64.9 60.0 5.01 5.5 

Experimental 6.241 — 



Other 6.34 b — 68 b 



42 


8.9 


37 


9.4 


35-38 c 


9.0 c 


35 b 


8.6 b 


21 


8.6 


20 c 


9.0 c 


20.8 b 


6.9 b 



a Ref.[7] 
b Ref. [5] 
c Ref [19] 
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Figure Caption 

Variation of total energy of PuP(a), PuAs(b) and PuSb(c) with relative volume in ferromagnetic 
and non-magnetic states in Bi phase. 

Variation of total energy of PuP(a), PuAs(b) and PuSb(c) with relative volume in Bi and B 2 
phases. 

variation of enthalpies in Bi and B 2 structures, and Pressure volume relation for PuP(a), 
PuAs(b) and PuSb(c). 

Band structure and Total DOS of ferromagnetic PuP(a), PuAs(b) and PuSb(c) at ambient 
pressure. Solid lines represent majority spin, doted lines represent minority spins. 

Band structure and Total DOS of ferromagnetic PuP(a), PuAs(b) and PuSb(c) at high 
pressure. Solid lines represent majority spin, doted lines represent minority spins. 
The partial density of states for PuP at the ambient pressure; Solid lines represent majority spin, 
doted lines represent minority spins. 

The partial density of states for PuAs at the ambient pressure; Solid lines represent majority spin, 
doted lines represent minority spins. 

Fig 8. The partial density of states for PuSb at the ambient pressure; Solid lines represent majority spin, 

doted lines represent minority spins. 
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Figure 3. 
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Figure 4 
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Figure 5 
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Figure 6 
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Figure 7 
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Figure 8 
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